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Abstract 

We analyze the energy spectrum of reactor neutrino events recently observed in the Kamioka 
Liquid scintillator Anti-Neutrino Detector (KamLAND) and combine them with solar and terres- 
trial neutrino data, in the context of two- and three-family active neutrino oscillations. In the 
1v case, we find that the solution to the solar neutrino problem at large mixing angle (LMA) is 
basically split into two sub-regions, that we denote as LMA-I and LMA-II. The LMA-I solution, 
characterized by lower values of the squared neutrino mass gap, is favored by the global data fit. 
This picture is not significantly modified in the ?>v mixing case. A brief discussion is given about 
the discrimination of the LMA-I and LMA-II solutions with future KamLAND data. In both the 
2u and 3u cases, we present a detailed analysis of the post-KamLAND bounds on the oscillation 
parameters. 

PACS numbers: 26.65.+t, 13.15.+g, 14.60.Pq, 91.35.-x 
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I. INTRODUCTION 



The year 2002 is likely to be remembered as the annus mirabilis of solar neutrino physics. 
On April 20, direct and highly significant evidence for v e flavor change into active states 
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On October 8, the role of solar neutrino physics in shaping modern 
science was recognized through the Nobel Prize jointly awarded to Raymond Davis, Jr., and 
Masatoshi Koshiba, for their pioneering contributions to the detection of cosmic neutrinos 
Finally, on December 6, clear "terrestrial" evidence for the oscillation solution to the 



solar neutrino deficit was reported by the Kamioka Liquid scintillator AntiNeutrino Detector 
(KamLAND), through the observation of long-baseline reactor 17 disappearanc e (l5| . The 
seminal idea of studying lepton physics by detecting solar 

0,1 

Il7l ] and reactor |l6l Il7l HE] 

neutrinos keeps thus bearing fruits after more than 50 years. 

The KamLAND observation of V P disappearance |l5j ] confirms the current interpretation 



of solar neutrino data Ifll. Ill . Iiol . EE 

neutrino mass and mixing [22J, |23[ , and 



2l| in terms of v e — > v^ r oscillations induced by 
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within the so-called large mixing angle (LMA 
solar neutrino data 24] , matter effects [IE HE in adiabatic regime 
dominate the dynamics of flavor transitions in the Sun (see, e.g 
spectral data appear to exclude some significant 
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In this paper we analyze the first KamLAND spectral data [15j and combine them with 

IE IE IE 0, IE EE EH , assuming two- and three-flavor oscillations 
in order to determine the surviving sub-regions of the LMA solution. 
In the analysis we include the CHOOZ reactor data [38j] and, in the 3z/ case, also the relevant 
constraints on the larger mass gap Am 2 3> 5m 2 , coming from the Super-Kamiokande (SK) 
atmospheric 3^| and KEK-to-Kamioka (K2K) accelerator (IE 41] neutrino experiments, 



, (K2] 
[21]- 1 



according to the approach developed in 

In the 2u case, we find that the inclusion of the KamLAND spectrum basically splits the 
LMA solution into two sub- regions at "small" and "large" 5m 2 , which we call LMA-I and 
LMA-II, respectively (the LMA-I solution being preferred by the data). Such regions are 
only slightly modified in the presence of 3z/ mixing, namely, for nonzero values of the mixing 
angle #13. We also present updated bounds in the 3z/ parameter space (5m 2 , 612, #13). 

The structure of this paper is as follows. In Sec. II we describe our analysis of the 
KamLAND data. In Sec. Ill we combine such measurements with solar and CHOOZ data, 
assuming 2v mixing. In Sec. IV we extend the analysis to 3z/ mixing, including the SK+K2K 
terrestrial neutrino constraints. We draw our conclusions in Sec. V. 



1 Our notation for the squared neutrino mass spectrum is (m\,m\,m\) = (— 5m 2 /2, +5m 2 , ±Am 2 ) as in 
[i^l, the sign of Am 2 being associated to normal (+) or inverted (— ) hierarchy. Unless otherwise noticed, 
normal hierarchy is assumed. The mixing angles 9ij are defined as in the quark case |^1E3- 
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II. KAMLAND DATA INPUT 



In our KamLAND data analysis, we use the absolute spectrum of events reported in |l5j , 
taken above a background-safe analysis threshold of 2.6 MeV in visible energy E? The 
events below such threshold might contain a significant component of geological z7 e ' s |43|, 
whose large normalization uncertainties are poorly constrained at present by the KamLAND 
data themselves [l5(- Above 2.6 MeV, a total of 54 events is found (including at most one 
possible background candidate), against 86.8 events expected from reactors [151 ] . 



The observed ener gy s pectrum of events is analyzed as in |21[ , with the following improve- 
ments. According to [l5j, we adopt: (1) Relative fuel components 235 TJ : 238 u ■ 239 p u : 241 Pu 
= 0.568 : 0.078 : 0.297 : 0.057; (2) Absolute normalization of 86.8 events for no oscillations; 
(3) Energy resolution width equal to 7.5%/ E 1 ^ 2 ; (4) Thirteen bins in visible energy above 
2.6 MeV, with 0.425 MeV width. The experimental and theoretical (oscillated) number of 
events in each bin are denoted as N^ xpt and N° sc , respecti vely (i = 1, . . . , 13). 

Concerning the systematic errors, the information in |l5j] does not allow to trace the 
magnitude and the bin-to- bin correlations of each component. We have then approximately 
grouped such uncertainties into two main components, labelled by the index k = 1, 2: 
(1) energy scale error; and (2) overall normalization error. The first (second) uncertainty 
basically shifts the spectrum along the energy (event) coordinate. More precisely, the effect of 
the energy scale error (k = 1) is estimated by taking the (symmetrized) fractional differences 
c\ between the values of iV° sc calculated with and without shifts of ±1.91% [l5[ in the 
true visible energy, for each point in the oscillation parameter space. The normalization 
uncertainty (k = 2) is here assumed to collect all the other error components listed in 
(lif . up to a total of about 6%: cf = 0.06. The systematic shifts c\ of the A° sc 's are 
then implemented through the "pull approach" described in [2£j, namely, through linear 
deviations N° sc -> A* heo = A° sc (l + £ fc c£&), where the &s are univariate Gaussian random 
variables. In this way, systematic correlations are also accounted for [201 ]. 

Concerning the statistical errors, the presence of bins with few or zero events requires 
Poisson statistics, that we approximately implement through the x 2 -like recipe suggested in 
the Review of Particle Properti es 42l | . A quadratic penalty in the pulls is introduced to 

The final x 2 function for KamLAND (KL) is thus 



account for the systematics |20j, |44 



Xkl = min 

t?fe> 



13 / yyexp \ 

j22(N^°-Nr+Nrin^)+ e e k 



fc=l,2 



(1) 



where the z-th logarithmic term is dropped if A 2 exp = |42j. By expanding the first sum in 
the x 2 a t first order in the shifts c^, the minimization in Eq. {1} becomes elementary. 

Some final remarks are in order. A detailed comparison of the above function Xkl with 
the one adopted by the KamLAND collaboration is not currently possible, since the latter 
is given in a symbolic form in 0|. However, one can observe that: (1) The KamLAND 
collaboration splits the total rate and spectrum shape information, while we prefer to use the 
absolute spectrum; (2) The KamLAND collaboration can take into account a (very small) 
background component above 2.6 MeV and four sources of systematics while, for a lack of 
information, we neglect such small background, and group the systematics into two main 
sources. For such reasons, we do not expect perfect quantitative agreement between the 



2 The "visible" or "prompt" energy E is denned as E = E v - (M n - M p ) + m e ~E u - 0.8 MeV. 
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TABLE I: Results of the two-flavor oscillation analysis. The first two columns report the data 
set used and the corresponding number of data point iVdata, respectively. The third and fourth 
columns give the position and x min value of the first and second best fit point (when applicable). 
For a goodness-of-fit check, the Xmin value should be compared with its expected ±lcr range, given 
by [A DF - y/2Nm,N DF + V^f] 0, where N DF = A data - 2 (see last column). All the x min 



values appear to be close to — or slightly below — the lower end of such ±lcr range, indicating that 
the fits are somewhat better than statistically expected (although not suspiciously so). 



Experimental data set 


No. of data 


Best fit point (s) a 


X 2 - 

A mm 


N DF ± V2A r DF 


Solar+CHOOZ 


81+14 


(5.5, 0.305) 


78.8 


[79.4, 106.6] 


KamLAND 


13 


(7.3, 0.335) 


6.1 


[6.3, 15.7] 




13 


(18.0, 0.270) 


7.9 


[6.3, 15.7] 


Solar+CHOOZ+KamLAND 


81+14+13 


(7.3, 0.315) 6 


85.2 


[91.4, 120.6] 




81+14+13 


(15.4, 0.300) c 


90.6 


[91.4, 120.6] 



"Coordinates are (Sm 2 /10" 5 eV 2 , sm 2 12 ). For KamLAND only, sin 2 6*i2 and 1 — sin 2 9\2 are equivalent. 
b Global best fit (LMA-I). 
c Second global best fit (LMA-II). 



KamLAND official oscillation analysis and ours. In fact, we find that the 95% C.L. contours 
of the rate+shape analysis in Fig. 6 of [lH are accurately reproduced by our Eq. (JTJ) at a 
slightly lower C.L. (90%). This loss of statistical power appears tolerable to us, in the light 
of the limited experimental information which is currently available. 3 



III. 2u ANALYSIS 

The updated 2z/ analysis of current solar+CHOOZ neutrino data, as performed in |2lj], is 
resented here for the sake of completeness. The fit includes 81 solar neutrino observables 



pres 
[20, 



20, 21 



and 14 CHOOZ spectrum bins |2lL 38j | . for a total of 95 data points. The best-fit 
point and its Xmin value are given in the first row of Table I. The A% 2 expansion around the 
minimum, relevant for the estimation of the oscillation parameters (5m 2 , sin 2 # 12 ), is shown 
in Fig. 1, where we have restricted the Sm 2 range to the only three decades relevant for the 
LMA solution and for the following KamLAND analysis. Notice that the scale is linear in 
the sin 2 9u variable. 

The 2z/ analysis of KamLAND is performed by using Eq. (0), which gives the absolute 
Xmin value reported in the second row of Table I. For later purposes, we also quote the second 
best fit parameters. Expansion around the absolute minimum gives the C.L. contours shown 
in Fig. 2. 4 There appears to be a "tower" of solutions which tend to merge and become 
indistinguishable for increasing Sm 2 ; the lower three ones are, however, rather well separated 
at 90% C.L. Notice that our allowed regions are slightly larger (i.e., less constraining) than 
those in the rate+shape analysis of [l5j, as explained in Sec. II. One of the two octant- 



3 A better evaluation and decomposition of systematic effects and an improved Xkl definition will be 
possible when more detailed KamLAND information on spectrum shape errors will become public. 

4 No oscillations (basically, oscillations with 5m 2 < 10~ 6 eV 2 ) provide a very bad fit, Ay 2 = 14.1. 
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symmetric best fits points in Fig. 2 (black dots) is remarkably close to the best fit in Fig. 1 
(see also Table I). The difference in location with respect to the KamLAND official best-fit 
point at (<5m 2 /eV 2 , sin 2 6*12) = (6.9 x 1CT 5 , 0.5) [l5| is not statistically significant, amounting 
to a variation A% 2 <C 1. 

The combination of the solar+CHOOZ results in Fig. 1 with the KamLAND results in 
Fig. 2 gives the global 2u results shown in Fig. 3, which represent the main result of this 
work. Two rather distinct solutions, that we label LMA-I and LMA-II, are seen to emerge. 
They are basically located at the intersection of the LMA solution in Fig. 1 with two of 
the well-separated KamLAND solutions in Fig. 2, and are characterized by the mass-mixing 
parameters and \ 2 values given in the last two rows of Table 1. The LMA-I solution is 
clearly preferred by the data, being close to the best fit points of both solar+CHOOZ and 
KamLAND data. The LMA-II solution is located at a 5m 2 value about twice as large as for 
the LMA-I, but is separated from the latter by a modest A% 2 = 5.4 difference (dominated 
by solar neutrino data). Indeed, if we conservatively demand a 99.73% CL. for the allowed 
regions, the LMA-I and LMA-II solutions appear to merge (and extend towards 5m 2 ~ 
3 x 10~ 4 eV 2 ) in a single broad solution. In any case, at any chosen C.L., the allowed regions 
of Fig. 3 are significantly smaller than those in Fig. 1. Therefore, with just 54 initial events, 
the KamLAND experiment is not only able to select the LMA region as the solution to 
the solar neutrino problem, but can also significantly restrict the corresponding oscillation 
parameter space. With several hundred events expected in the forthcoming years, there are 



thus very good prospects to refine the (5m 2 , sin 2 # 12 ) parameter estimate [21 

The most important task for the next future appears to be the confirmation of one of 
the two solutions in Fig. 3 through higher statistics and, possibly, lower analysis threshold. 
Figure 4 shows the absolute KamLAND energy spectra of reactor events predicted at the 
LMA-I and LMA-II global best-fit points (last two rows of Table I), together with the no- 
oscillation spectrum, with the current event normalization. The KamLAND data (with 
statistical errors only) are also superposed above the analysis threshold (2.6 MeV). It can 
be seen that the main difference between the two oscillated spectra is the position of the 
spectrum peak, which is roughly aligned with the no-oscillation position for the LMA-II 
case, while it is shifted at higher energies for the LMA-I case. This feature might be a useful 
discrimination tool in the next future. From Fig. 4, it appears also that the largest difference 
between the LMA-I and LMA-II spectra in KamLAND occurs just in the first bin below 
the current analysis threshold. Therefore, a better understanding of the background below 
2.6 MeV will be highly beneficial. In any case, an increase in statistics by a factor of a few 
appears necessary to disentangle the two oscillated spectra in Fig. 4. 

IV. 3i/ ANALYSIS 

The updated 3i> analysis of solar+CHOOZ neutrino data, including the constraints on 
the large ("atmospheric") squared mass gap Am 2 coming from SK atmospheric j3^| and 
K2K accelerator neutrino data ji^, 41|, as performed in j2lj], is reported here for the sake 
of completeness. 5 Figure 5 shows the pre-KamLAND 3u results 21] in the parameter space 



In [21| , the combination of SK and K2K constraints on the Am 2 parameter was obtained by adding the 
preliminary x 2 functions from the two collaborations, as graphically presented in |sj],^(J. We anticipate 
that such combination appears to be in good agreement with a more refined and joint reanalysis of SK 
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(5m 2 , sin 2 #12, sin 2 $13) relevant for the LMA solution, shown through sections at four different 
values of s 2 3 = sin 2 #13. 6 For later purposes, we only notice that the upper bound on 5m 2 
becomes stronger for higher sin #13, as a result of the CHOOZ constraints. For a discussion 
of such (Sm 2 , sin 2 # 13 ) anti-correlation, and for the pre-KamLAND solutions below the LMA, 
see |2lJ and references therein. 

Figure 6 shows the results of our 2>v fit to the KamLAND data, in the same format as 
in Fig. 5. The allowed 3z/ regions appear to get slightly en larg ed (especially in the mixing 
parameter #i 2 ) for increasing values of sin 2 #13, as expected |34j. In fact, for #13 7^ 0, part of 
the KamLAND event disappearance is explained by averaged oscillations driven by Am 2 in 
the (v e , i/ 3 ) sector, so that the overall oscillation amplitude sin 2 2#i 2 in the (^1,^2) sector is 
allowed to reach smaller values, and the sin 2 Qyi range is correspondingly enlarged [2H . The 
absolute x 2 minimum in Fig . 6 is reached for 6* 13 7^ 0, but such preference is not statistically 
significant, as expected [30(; indeed, we find a mere A% 2 = 0.4 increase of the KamLAND 
Xmin vame f° r fixed sin 2 9 13 = 0.06. Concerning the variations of the (5m 2 , sin 2 # 12 ) best-fit 
coordinates for fixed values of sin 2 #13, we find that the Sm 2 value is stable (and equal to 
7.3 x 10~ 5 eV 2 , as in the 2z/ case), while sin 2 # 12 decreases from 0.335 to 0.225 when increa sing 
sin 2 6*13 from to 0.06, in good qualitative agreement with the expectations discussed in J34J. 

However, the additional spread in sin 2 # i2 induced by nonzero 6*13 currently does not 
play any relevant role when KamLAND is combined with world neutrino data, for at least 
two reasons: (1) Pre-KamLAND bounds from solar+terrestrial data currently dominate the 
constraints on sin 2 # 12 , as evident from a comparison of Fig. 5 and Fig. 6; (2) The likelihood 
of genuine 3i/-induced effects rapidly decreases with increasing 13 , because of the strong 
upper bounds on such mixing parameter [2l|. Therefore, we do not expect any significant 
enlargement of the sin 2 # 12 allowed range from the pre- to the post-KamLAND 3u analysis, 
despite the presence of such effects in KamLAND alone. 

It should be noted that, strictly speaking, Fig. 6 is not exactly a "KamLAND only" 
analysis, since we have implicitly taken some pieces of information from terrestrial neutrino 
data. In particular, we have implicitly assumed in Fig. 6 that: (i) the "atmospheric" squared 
mass splitting Am 2 is sufficiently high to be unresolved in KamLAND, and (2) the relevant 
values of sin 2 #13 are limited in the few percent range. It is interesting to study the effect 
of an explicit combination of such (atmospheric + CHOOZ) information with KamLAND 
data. The results are given in Fig. 7. This figure show that purely terrestrial neutrino data 
from atmospheric (SK), accelerator (K2K) and reactor (KamLAND + CHOOZ) neutrino 
experiments, by themselves, are now able to put both upper and lower bounds on the solar 
parameters (5m 2 , sin 2 812). As previously noted in the comment to Fig. 5, the CHOOZ upper 
bound on 5m 2 becomes stronger when sin 2 #i 3 increases. 7 The slight octant asymmetry for 
sin 2 #i3 > is due to the fact the corresponding 3u CHOOZ probability is not invariant under 
the change #i 2 — > ir/2 — # 12 for fixed hierarchy (assumed to be normal in Fig. 6). The 
octant differences would be swapped by inverting the hierarchy [471] . In practice, however, 
such octant asymmetries are numerically irrelevant in the global fit which we now discuss. 

Figure 8 shows the final 3z/ combination of world (solar+terrestrial) neutrino con- 
straints, including pre-KamLAND data (Fig. 5) and the first KamLAND data (Fig. 6). 



atmospheric and K2K data [4f 

6 We remind that, in Fig. 5 and in the following iv figures, the information on Am 2 is projected away by 
taking minA m 2(x 2 ), see [ilj . 

7 Analogously, the upper bound on sin 2 #13 becomes stronger for increasing 5m 2 |46( . 
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For sin 2 #13 = (absolute best fit), we get the same LMA-I and LMA-II solutions reported 
in Fig. 3 and in Table I, modulo the expected widening induced by one additional degree 
of freedom in the A% 2 's associated to each C.L. contour. This widening gives, at 99% C.L., 
marginal allowance for a third solution with 5m 2 ~ 2.5-3.2 x 10~ 4 eV 2 , which we call LMA- 
III. 8 All solutions rapidly shrink for increasing sin 2 # 13 , the LMA-I being the most stable 
and the last to disappear. The hierarchy is assumed to be normal in Fig. 8; the differences 
with respect to the inverted hierarchy case (not shown) are completely negligible. 

Figure 9 represents a compendium of our post-KamLAND 3u analysis of all data 
(solar+terrestrial), in terms of the A% 2 function for each of the three variables 
(5m 2 , sin 2 #i2, sin 2 #i 3 ), the others being projected (minimized) away. This figure should be 
compared with the pre-KamLAND one in |2l|. The comparison shows that the bounds on 
sin 2 #13 remain basically unaltered in the post-KamLAND era, while the mixing parameter 
sin 2 #12 is clearly nailed down around the value 0.3. The absolute and second best minima 
(LMA-I and LMA-II) are also evident in the left panel. A third minimum (LMA-III) at 
higher 5m 2 is only marginally allowed. 

If we focus on the best-fit (LMA-I) solution, and adapt parabolic functions around the 
absolute x 2 minima in terms of 5m 2 (in linear scale), sin 2 #13, and sin 2 #13, we obtain from 
the post-KamLAND global analysis the following approximate ±1<t estimates (Ax 2 ~ 1) for 
the relevant solar 3z/ oscillation parameters, 

f Sm 2 ~ (7.3 ± 0.8) x lO" 5 eV 2 , 
LMA-I (~ la) : J sin 2 # 12 ~ 0.315 ± 0.035 , (2) 
[ sin 2 #13 < 0.017 . 

The above ranges are meant to show that we are not far from a 10% determination (at 
la) of both 5m 2 and sin 2 #i 2 , but should not be quoted as a "summary" of the current 3u 
situation. The technically correct reference summary for our analysis is represented by the 
X 2 functions of Fig. 9, which include the possibility of a second global best fit (LMA-II), 
and maybe of a third best fit (LMA-III) in Sm 2 . In particular, we stress that the LMA-II 
solution is currently quite acceptable from a statistical point of view, and that the history of 
the solar neutrino problem teaches us that one should not exclude a priori what may appear 
as the second best-fit at a particular time. 

Similar cautionary remarks also apply to the use of the initial KamLAND data to refine 
the emerging indications of adiabatic matter effects in the Sun 2^, 2?} through the approach 



developed in [28[. As stressed in [28[, such analysis will become sufficiently stable and 
constraining when a single solution (either the LMA-I or the LMA-II) will be (hopefully) 
clearly selected by KamLAND. In the meantime, further SNO data are also expected to 
increase our confidence in the occurrence of solar matter effects, see |28| . 

We conclude our 3z/ analysis by observing that, from the point of view of future (very) 
long baseline (LBL) experiments, the "solar" neutrino parameters 5m 2 and #12 enter the 
subleading oscillation amplitudes of both the "golden" channel v e — > j48| and the "silver" 
channel v P — > v T ji^] only through the combination n = 5m 2 sin2#i2, up to 0(n 2 ) terms in- 
cluded |48l.l49l]. Therefore, we think it useful to show in Fig. 10 the global (solar+terrestrial) 
bounds on this specific combination of parameters (in units of 10 -5 eV 2 and in linear scale), 
all other 3u variables being projected away. The lower and higher minima in Fig. 10 cor- 
respond to the LMA-I and LMA-II cases, respectively. By adapting a parabola in the ±2<r 



Notice that the LMA-I and LMA-III solutions differ by Ax 2 ^ 11. 
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region around each of the two minima, we get the following la (A% 2 = 1) estimates for 
re = 5m 2 sin 29 12 : re ^ 6.7 ± 0.7 (LMA-I) and re ~ 14.7 ± 1.7 (LMA-II), in units of KT 5 
eV 2 . In both cases, re appears to be already constrained with a ~ 11% accuracy. These 
approximate estimates — or, better, the accurate bounds in Fig. 10 — may prove useful in 
prospective studies of re-dependent subleading oscillation effects in LBL projects. 

V. SUMMARY AND PROSPECTS 

The KamLAND experiment has clearly selected the LM A region as the solution to the 
solar neutrino problem, and has further reduced the (5m 2 , sin 2 9 12 ) parameter space for active 
neutrino oscillations. In this work, we have studied such parameter reduction in the the 
context of 2is and 3z/ oscillations, by using the limited KamLAND experimental information 
which is publicly available. In the 2is case, we find that the post-KamLAND LMA solution 
appears to be basically split into two sub-regions, LMA-I and LMA-II. The LMA-I solution, 
characterized by 5m 2 ~ 7 x 10 -5 eV 2 and sin 2 9 12 ~ 0.3, is preferred by the global fit. 
The LMA-II solution represents the second best fit, at about twice the value of 5m 2 (see 
Table I). This situation is not significantly changed in the 3z/ case, for which we present a 
global post-KamLAND analysis of solar and terrestrial data in the (5m 2 , sin 2 # 12 , sin 2 9 13 ) 
parameter space. There are good prospects to separate the LMA-I and LMA-II cases with 
future, higher-statistics KamLAND data, by looking at the peak of the energy spectrum, 
and by lowering the current analysis threshold (2.6 MeV) by at least ~ 0.5 MeV. 
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FIG. 1: Two-flavor active neutrino oscillations: Global analysis of solar and CHOOZ neutrino 
data in the (5m 2 , sin 2 612) parameter space, restricted to the LMA region. The best fit is indicated 
by a black dot. The 90, 95, 99, and 99.73% C.L. contours correspond to A X 2 = 4.61, 5.99, 9.21, 
and 11.83, respectively (two degrees of freedom). For solutions at lower 5m 2 — now ruled out by 
KamLAND— see E3. 
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FIG. 2: Two- flavor active neutrino oscillations: Analysis of the KamLAND energy spectrum data 
above 2.6 MeV in the (5m 2 , sin 2 612) parameter space. A "tower" of octant-symmetric regions is 
allowed at different values of 5m 2 . The three lower regions with 5m? < 2 x 10 -4 eV 2 are well 
separated at 90% C.L., while the upper ones tend to merge in a continuum. The symmetric best 
fits are indicated by black dots. The left dot is remarkably close to the solar best fit in Fig. 1. 
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2v Solar + CHOOZ + KamLAND 




FIG. 3: Two-flavor active neutrino oscillations: Global analysis of solar, CHOOZ, and KamLAND 
neutrino data in the (5m 2 , sin 2 612) parameter space. With respect to Fig. 1, the LMA region is 
significantly restricted, and appears to be split into two sub-regions (LMA-I and LMA-II), well- 
separated at 99% C.L. The best fit (black dot) is found in the LMA-I solution. At 99.73% C.L., 
the two solutions merge into a single one, which slightly extends towards 5m 2 ~ 3 x 10~ 4 eV 2 . 
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FIG. 4: Energy spectra of events in KamLAND for the LMA-I and LMA-II global fit points in 
Table II, together with the no-oscillation spectrum. The current KamLAND data are superposed 
above the analysis threshold (2.6 MeV). 
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3p oscillations : Solar + (Atm. + K2K) + CHOOZ 
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FIG. 5: Three-flavor active neutrino oscillations: Global analysis of solar, CHOOZ, and atmo- 
spheric+K2K neutrino data, shown as sections of the allowed (5m 2 , sin 2 #12, sin 2 #13) parameter 
space at four representative values of s 2 3 = sin 2 #13. The best fit is reached in the 2v limit 
(s 2 3 = 0), while for increasing s 2 3 the LMA solution shrinks and eventually disappears. The 90, 
95, 99, and 99.73% C.L. contours correspond to A% 2 = 6.25, 7.82, 11.34, and 14.15, respectively 
(three degrees of freedom). The "atmospheric" mass splitting Am 2 is minimized away in the fit 
(see also |2l|). 
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FIG. 6: Three-flavor active neutrino oscillations: Global analysis of KamLAND data, shown in 
the same format as Fig. 5. Notice the enlargement of the allowed regions for increasing sf 3 . 
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FIG. 7: Three- flavor active neutrino oscillations: Global analysis of terrestrial data only (including 
KamLAND), shown in the same format as Fig. 5. Notice that terrestrial v data alone can now put 
both upper and lower bounds on the solar v parameters (5m 2 , sin 2 612)- See the text for details. 
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3u oscillations : Solar + (Atm. + K2K) + CHOOZ + KamLAND 
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FIG. 8: Three-flavor active neutrino oscillations: Global analysis of solar and terrestrial data 
(including KamLAND), shown in the same format as Fig. 5. With respect to the 1v case in Fig. 3, 
the higher Ax 2 tolerance induced by an extra degree of freedom (for any chosen C.L.) marginally 
allows a third "LMA-III" solution at 99% C.L., at about Sm 2 ~ 2.5-3.2 eV 2 . 
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FIG. 9: Three-flavor active neutrino oscillations: Projections of the global (Ax 2 ) function onto 
each of the (dm 2 , sin 2 9±2, sin 2 #13) parameters. The n-sigma bounds on each parameter (the others 
being unconstrained) correspond to A% 2 = n 2 . In the left panel, the absolute minimum (LMA-I) 
and the second best fit (LMA-II) are clearly visible. A third shallow minimum (LMA-III) of the 
Ax 2 (Sm 2 ) function is marginally present at ~ 3.2a. 
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FIG. 10: Global bounds on the parameter combination 5m 2 sin 28±2 , governing the amplitude 
of subleading "solar" neutrino oscillations in long baseline (LBL) experiments. See the text for 
details. 
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